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Abstract 

By the introduction of a generalized Evans function defined by an appropriate 
2-modified Fredholm determinant, we give a simple proof of convergence in location 
and multiplicity of Hill's method for numerical approximation of spectra of periodic- 
coefficient ordinary differential operators. Our results apply to operators of nondegener- 
ate type, under the condition that the principal coefficient matrix be symmetric positive 
definite (automatically satisfied in the scalar case). Notably, this includes a large class 
of nonselfadjoint operators, which were previously not treated. The case of general 
coefficients depends on an interesting operator-theoretic question regarding properties 
of Toeplitz matrices. 

1 Introduction 

The study of stability of spatially periodic traveling wave solutions to various classes of 
partial differential equations motivates the study of L 2 (M; C n ) (essential) spectra of periodic- 
coefficient differential operators 

(1.1) L = (d x ) m a m (x) + ■■■ + d xai {x) + a (x) 

on the line, where coefficients aj £ C raxn are periodic with period X. By Floquet theory, it 
is equivalent to study the L 2 ([0, X];C n ) point spectra of the family of Bloch operators 

L a = (d x + ia) m a m (x) H h (d x + ia)ai{x) +a (x), 
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where X is the common period of the coefficients and a G [0, 2tt) acts as a parameter. 
Indeed, using this decomposition we hav^] 



see, for example, [G] for more details. 

Due to the mathematical difficulties involved in analytically computing the L 2 (R) spec- 
trum of such an, in general, variable-coefficient and vector- valued, operator, or, equivalently, 
computing the periodic spectra of the full family of associated Bloch operators, the deter- 
mination of spectrum of periodic-coefficient operators is typically carried out numerically. 
This may be accomplished in a number of ways: for example, shooting, discretization, or 
various spectral and Galerkin methods. See Appendix B, [JZNj . for further discussion. 



method carried out in a periodic Fourier basis, which is exact in the constant-coefficient 
case. In this method, to approximate the spectra of L a for a fixed a € [0, 2tt), one considers 
the eigenvalue problem 



by expressing the coefficients ctj of L a and the function v as Fourier series in Lp er ([0, X]), as 
an infinite-dimensional matrix equation in £ 2 . Truncating the Fourier modes to frequencies 
\k\ < J for each J £ N, one then obtains a sequence of finite-dimensional matrix eigenvalue 
problem whose eigenvalues approximate true eigenvalues of the operator L a on Lp er ([0, X]). 
See Section [3T21 for further details. 

This method is fast and easy to use, and in practice appears to give excellent results 
under quite general circumstances [DK} iBJNRZl], However, up to now, an accompanying 
rigorous convergence theory has been established only in certain commonly occurring but 
restricted cases |CuD| . By convergence, we mean roughly that not only is Hill's method 
accurate, meaning that the numerically computed eigenvalues are always close to the actual 
eigenvalues of the associated Bloch-operator (the "no-spurious modes condition" of [CuDj ) , 
but also that the method is complete in the sense that it faithfully produces all of cr(L a ) for a 
fixed a: see [CuD for a more precise discussion of convergence from this point of view. Here, 
we make the simpler, operational definition that on any bounded domain B = {A : |A| < R} 
whose boundary contains no eigenvalue of L a , the set of approximate eigenvalues lying in 
B converges to the set of exact eigenvalues of L in both location and number; see Cor. 13.9^1 

Despite its obvious practical interest, up to now the convergence of Hill's method has 
been established to our knowledge only for self-adjoint operators with principal coefficient 
a m = I [CuD] , In particular, though accuracy of Hill's method was shown in |CuD| under 
quite general assumptions, completeness of the method in the non-selfadjoint case, which 

1 Unless otherwise stated, throughout this paper all functions are assumed to be complex valued and we 
adopt the notation L 2 (R) = L 2 (R;C) and similarly for Lp er ([0, X]). 

2 A convenient implementation may be found in the numerical package SpectrUW [CDKK| . 
3 This includes and slightly strengthens the definition of [CuD] . 
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A particularly natural and direct approach is Hill's method 




(1.2) 



L a v = Xv 
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arises naturally, for example, in the applications in [BJNRZ1, BJNRZ2], does not seem to 
have been fully addressed. 

In this short paper, we give a brief and simple proof of the convergence of Hill's method 
applying to the general class of operators (jl.ip such that a m is symmetric positive defi- 
nite. In the scalar case, this condition on the principal coefficient a m amounts to the mild 
requirement that the operator be nondegenerate type. In the system case, it is a genuine 
restriction, and it is an interesting and apparently nontrivial question, related to certain 
properties of Toeplitz matrices, to what extent the condition can be relaxed. Notably, our 
analysis applies to the important case where the operator L a is non-selfadjoint. 

The main ingredient of our our proof is the introduction of a generalized periodic Evans 
function, of interest in its own right, consisting of a 2-modified Fredholm determinant D a 
of an associated Birman-Schwinger type operator, whose roots we show to agree in location 
and multiplicity with the eigenvalues of L a . For related analysis in the solitary wave case, see 
[GLZ]. Once these properties are established, the desired convergence follows immediately 
by the observation that the corresponding 2-modified characteristic polynomial of the J th 
Galerkin-truncation of (L a — X)v = are a subclass of the approximants used to define the 
aforementioned 2-modified Fredholm determinant in the limit as J — > oo, and furthermore 
that these approximates are a sequence of analytic functions converging locally uniformly 
to the generalized periodic Evans function. 

A novel feature of the present analysis is that our argument yields convergence of the 
spectrum in both location and multiplicity, whereas the results of [CuD] concerned only 
location. On the other hand, there was established in [CuD] a fast rate of convergence to 
the smallest (in modulus) eigenvalue in the self-adjoint case, whereas our methods do not 
readily appear to yield a rate. A second novelty of our work is to make the connection to 
the Evans function, putting the work in a broader context. 

2 Hilbert— Schmidt operators and 2-modified Fredholm de- 
terminants 

We begin by recalling the basic properties of 2-modified Fredholm determinants, defined 
for Hilbert-Schmidt perturbations of the identity; see [GGKH IGGK2] . [GGK31 Ch. XIII], 
[GKl Sect. IV.2], fSH], [SHI Ch. 3] [GLZl Sect. 2] for more details. 

For a given Hilbert space %0 the Hilbert-Schmidt class B2CH) is defined as the set of 
all bounded linear operators A on % for which the norm 

\\ A \Wu) -=Y.\( Ae ^ ek )\ 2 = tr «( AM ) 

is finite, where {ej} is any orthonormal basis. Evidently, || • ||b 2 (W) i s independent of the 
basis chosen. Moreover, every operator in B2CH) is compact (Fredholm). 



4 Throughout this paper, we will always assume that our Hilbert spaces are separable. 
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On a finite-dimensional space H, we define the 2-modified Fredholm determinant as 
(2-1) det 2iW (J« - A) := det H ((I n - A)e A ) = det n (I n - A) e tr ^ A \ 

where det% and tr% denotes the usual determinant and trace, respectively. From this 
definition, we have the useful estimates 

(2.2) |det 2 ,«(/« - A)\ < e C!|Ai| s2(«) 
and 

(2.3) |det 2 ,«(/« - A) - det 2 , H ( I H ~ B)\ < \\A - S|| e2(w)e c; [!l A !l«2(^)+ll B ll82(«)+ 1 ] 2 , 

where C > is a constant independent of the dimension of %. 

To extend this notion of a determinant to an infinite dimensional Hilbert space T~L, we 
note that for any A G £> 2 (%) the estimate (j2.3|) allows us to define the 2-modified Fredholm 
determinant unambiguously as the limit 

(2.4) det 2 ,n{In ~ A) := lim det 2 , Wj (/ W/ - Aj), 

J— ¥00 

where %j is any increasing sequence of finite-dimensional subspaces filling up H, and Aj 
denotes the Galerkin approximation Pu A\nj, where Pj : % — > T-Lj is the orthogonal 
projection onto T~Lj. That is, thinking of the infinite-dimensional matrix representation of 
A, the 2-modified Fredholm determinant is defined as the limit of such determinants on 
finite, J-dimensional, minors as J — > oo. 

Alternatively, denoting the (countably many, since A is Fredholm) eigenvalues of A as 
{aj}'jL l , and taking %j to be the (total) eigenspace associated with the eigenvalues {aj}j =1 
we find that 

J 

(2.5) det 2 ,w(/w — A) = lim TT(1 - a k )e a \ 

J— >oo 

fc=l 

which, by Hfc(l — ctk)e ak < 11^(1 + a\) ~ e^ fcQ fe < e"" 4 "^^), is readily seen to converge 
for all A G ;6 2 (%) by Weyl's inequality ^l a j| r — Sl s j| r for r > 0, where Sj denote 
the eigenvalues of \A\ := (A* A) 1 / 2 jSil|. IW| . This shows how the renormalization of the 
standard determinant det(/% — ^4) := IIj(l — ay) by factor e^ lH ^ cancels the possibly 
divergent first-order terms in 14^(1 — a k ) ~ e^-- feQfe , allowing the treatment of operators A 
that are not in trace class B\ := {A : IWA^^W^^ < +oo}[j 

Proposition 2.1. For A G £> 2 (%), the operator (/% — A) is invertible if and only if 
det 2i %(/% — A) is non-zero. 

5 For A £ Bi, tr-H(A) — £\ ctj is absolutely convergent, by Weyl's inequality with r = 1, and so the 
standard determinant det-«(/^ — A) = Hj(l — aj) converges. For A self-adjoint, II^Hbj := || l-AI 1 ' 2 ||b 2 (-h) = 
Ej \ a j\ and \\ a \\b 2 (H) = Ej \ a A 2 - 
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Proof. By standard Fredholm theory, this is equivalent to the statement that is an eigen- 
value of (In — A) if and only if det2,u(I'H — ^4) = 0. Note that, since A is Fredholm, it 
possesses a countable number of isolated eigenvalues {<x/} of finite multiplicity, except pos- 
sibly at zero. Choosing J 6 N sufficiently large, then, we may factor the product formula 
flZSD as 

det 2 , n (lH-K)= I J ( (l-^)e C 

\i=i / \i=J+i 

where 

oo 

[J (1 - o^e"' ~ e E ^+i a ? / 0. 
j=J+i 

It follows then that det2 1 n(In — A) vanishes if and only if 1 — ay = for some 1 < j < J, 
hence, since J £ N was arbitrary, if and only if is an eigenvalue of (/% — A). □ 



3 Analysis of a simple case 

With the above preliminaries in hand, we now turn to our proof of convergence. As a first 
step in this analysis, we present a complete proof in the case of a second-order operator with 
identity principal part. In later sections, we will then describe the extension of this proof 
to more general cases, noting that most of the ideas can be found in this simpler context. 
Consider a periodic-coefficient differential operator 

La = (d x + io~) 2 + (d x + icr)ai(x) + a (x) 

acting on vector-valued functions in Lp er ([0,X]), a £ [0, 2tt) the Floquet parameter and 
cij G L 2 ([0, X]) matrix- valued and periodic on x £ [0, X]. We can rewrite this more generally 
as a family of operators in the simpler form 

(3.1) L a = dl + d x A 1 (a,x) + A (a,x), 
where 

A\ = ai + 2ia, Aq = (iq — a 2 + iaa\. 

In order to analyze the (necessarily discrete) spectrum of the operator L a , we introduce 
a generalization of the periodic Evans function, a complex analytic function whose roots 
coincide in location and multiplicity with the eigenvalues of L a [G] , expressed in terms of a 
2- modified Fredholm determinant. To this end, notice that associated with the eigenvalue 
problem 

(3.2) (L a -X)U = 
is the equivalent problem 



(3.3) 



(I + K(a,X))U = 0, 
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where here I is the identity operator on L2 er ([0, X}) and K = K x + K , with 

K x = d x {d 2 x - l)- l A x , K = (d 2 x - l)-\A + 1 - A). 

In particular, notice that A is an eigenvalue of L a if and only if is an eigenvalue of the 
operator (/ + K(a, A)). Before we can define the appropriate generalization of the Evans 
function, we need the following fundamental lemma. 

Lemma 3.1. For Aj € Lp CI .([0,X]) ; the operator K is Hilbert- Schmidt. 

Proof. Expressing K m in matrix form fC m with respect to the infinite-dimensional Fourier 
basis, we find that the corresponding matrix elements can be expressed as 

where A\(m) denotes the m th Fourier coefficient of A x , and i := \f—l. Computing explicitly, 
we find by Parseval's Theorem thatlfl 

;2 



\K x \\b 2 = ||/Ci||ft =Y, dl^ Yj^i 

j ^ J 1 k 
j2 

= S ( 1+J -2) 2 H A llU|e r (M < +°°> 



hence K\ is a Hilbert-Schmidt operator. Similarly, we find that Kq is Hilbert-Schmidt, 
with norm 

Uplift = E JTTpy S | io ^' " fc ) + ( x " f . 

j k 

which implies that K = K\ + Kq S B2 as claimed. □ 

Remark 3.2. On the other hand, K\ is not trace class if A±(0) := J A\(x)dx 7^ ; since 
then Y2j l^ljjl = 1-^1 (0)| Y2j = +°°- This illustrates the necessity of our extension of 

the usual notion of a determinant to operators in B2. 

3.1 Generalized Periodic Evans Function 

By Lemma 13. II in conjunction with Proposition 12. 1( it follows that the zero eigenvalues of 
(i£2 er ([o,X]) — K{a, A)) can be identified through the use of a 2-modified Fredholm determi- 
nant. This leads us to the following definition. 

Definition 3.3. For a fixed a £ [0, 2tt), we define the generalized periodic Evans function 
Dp '. (C — y C by 

(3-4) D a (X) := det 2iL 2 or([0iX]) (/ L 2 or([0iX]) - K(a, A)). 

6 Henceforth, Hilbert-Schmidt spaces B2 will always be considered on the Hilbert space Lp Cr ([0, X]). That 
is, we adopt the notation B2 ■= B2(Lp el ([0, X])). 
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For ease of notation, throughout the rest of our analysis we will drop the dependence on 
the Hilbert space Lp er ([0,X]) on the identity operator and all 2-modified Fredholm deter- 
minants. In particular, we will write D a (\) = det2(I — K(cr,X)) for the above generalized 
Evans function. 

Theorem 3.4. For Aj G L^ er ([0, X]), the function D a is complex- analytic in A and con- 
tinuous in the parameter a. Furthermore, the roots of D a for a fixed a G [0, 2tt) correspond 
in location and multiplicity with the eigenvalues of La. 

Proof. Following the notation in Lemma I3TTI for each J G N we let ICj := ([K]j,k)\j\,\k\<j be 
the finite dimensional Galerkin matrix approximation of the bi-infinite dimensional matrix 
representation of the operator K defined above. Clearly, then, for each fixed J G N the 
finite-dimensional approximation Aj(<t, A) := det2(/ — K.j(a,\)) is complex-analytic in A 
and continuous in a G [0, 2ir). Furthermore, as in the proof of Lemma 13.11 we have 



where JC\ t j denotes the truncation of /Ci, and hence we find that K\ t j — > K\ in B2 uniformly 
in both a and A. Similarly, we find that /Co,j(<7, A) — > /Cn(c, A) in B2 uniformly in a and 
locally uniformly in A, and hence the estimate f)2.3f) implieqll that Aj — > D a locally uniformly 
in A G C and uniformly in a G [0, 2tt). It follows that the function (a, A) 1— > D a (\) inherits 
the same regularity properties in A and a as the limiting sequence Aj, thus verifying the 
first claim of the Theorem. 

Next, by equivalence of the problems (j3.2j) and (j3.3[) together with Proposition 12. 1\ we 
obtain immediately correspondence in location of the roots of D a and the eigenvalues of 
the operator L a . To obtain agreement in multiplicity, consider an eigenvalue A* of L a , with 
corresponding eigenspace H*. Recalling that, by standard Fredholm theory, the eigenvalues 
of L a are countable, isolated, and have finite-multiplicit}{l, we find that there exists a closed 
ball -B(A*,e) of radius e, centered at A*, containing no other eigenvalues of L a . 

Consider now an increasing sequence of eigenspaces {-HjjjgN of Lp er ([0, X]) such that 
YvmjHj = Lp er ([0,X]) and C Hj for all J G N. For each J, let {r fc }^ =1 be an or- 
thonormal basis of Hj and let Rj = (r±, . . . ,rj). Then we can define the finite-dimensional 
approximants 



Since D a does not vanish on dB(\*,e), by the correspondence in location of roots and 
eigenvalues established above, and since Sj converges locally uniformly in A to D a by (|2.3|) . 
Rouche's Theorem implies that there exists a J* G N sufficiently large such that for J > J* 

7 To use the estimate (|2.3|) directly, one should consider the operator Kj, which is technically defined 
on the finite-dimensional subspace Hj, as being defined on the larger space Lp el ([0, X]). Throughout the 
remainder of our analysis we will consider this extension without reserve. 

8 Note that in this standard theory, one inverts L„ — [il rather than D 2 — 1. 



(3.5) 




lil>J+i 



(3.6) 
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the winding number of D a around dB(X*,e) is equal to the winding number of 5j around 
the same ball. 

Finally, fixing Jo > J* and noticing that L a Rj = Rj M a j , where M a j Q is an Jo x Jo 
matrix representation of L a on the finite-dimensional invariant subspace Hj , we find from 
(|3.6|) that there exists a constant C / such that 

S Jo (a, A) = det 2 {R* Jo (d 2 x - l) _1 i?j (M (Ti j - XI)) = Cdet 2 (M Jo - XI) , 

and hence we see that 8j is a nonvanishing multiple of the characteristic polynomial of 
M CTi j . Here, we are using the fact that Rj Q (d 2 — l)~ x Rj is positive definite, by positive 
symmetric definiteness of (d 2 — 1) . It follows that 5j has a zero at A* of precisely the 
algebraic multiplicity of A* as an eigenvalue of L a . Thus, we conclude that the multiplicity 
of A* as a root of D a is equal to the winding number of Sj (-,a) about the ball dB(X*,e), 
which in turn is equal to the algebraic multiplicity of A* as an eigenvalue of L a , completing 
the proof. □ 

Remark 3.5. The truncated winding-number argument for agreement of multiplicity to our 
knowledge is new, and seems of general use in similar situations. It would be interesting to 
prove this also in a different way by establishing a direct correspondence between the Fred- 
holm determinant and the standard periodic Evans function construction of Gardner 
as done in the solitary-wave case in \GLM1\ I GLM Z2. GMj and in the periodic Schrddinger 
case in [GM, Sect. 4]- This would give at the same time an alternative proof of Gardner's 
fundamental result of agreement in location and multiplicity of roots of the standard periodic 
Evans function with eigenvalues of L a , through the result of Theorem \3.4\ 

3.2 Convergence of Hill's method 

Next, we use the machinery developed in the previous section to give a proof of the conver- 
gence of Hill's method. In order to precisely describe Hill's method, notice that by taking 
the Fourier transform, we may express f|3.2[) equivalently as the infinite-dimensional matrix 
system 

(V 2 + VAi + Ao - XI)U = 0, 
where for each m = 0, 1 and j, k € Z, 

(3.7) V jk = 5 k j ij, [Am]jk = Anti ~ k), &ndU j = U(j), 

where f(k) denotes the discrete Fourier tranform of / evaluated at Fourier frequency k and, 
as elsewhere, i = \f— 1. Hill's method then consists of fixing J E N and truncating the above 
infinite-dimensional matrix system at wave number J, that is, considering the (2J + 1)- 
dimensional minor k)\ < J, and computing the eigenvalues of the finite-dimensional 
matrix 



(3.8) 



L CTiJ :=Vj + VjAi,j + A ,j, 
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where T>j and A m ^j denote the (2J + l)-dimensional matrices resulting from truncating 
the matrices V and A m to frequencies k)\ < J, to obtain approximate eigenvalues for 
L a . Notice this can be done quite efficiently by applying modern numerical linear algebra 
techniques. 

Remark 3.6. In applications, one may of course encounter operators L that are not in 
divergence form (]3.ip . In this case, we point out that there is no effect in changing from 
nondivergence to divergence form except that we increase the regularity requirement on A\ 
from L 2 to H 1 . Indeed,, we may change from one form to the other using the Leibnitz rule 
AiV - VA\ = (Ai)', where 

(Ai)' jk = i(j - k)At(j -k) = (A^ x )(j - k), 

and noting that, since T> is diagonal, this operation is respected by truncation. Thus, there 
is indeed no loss of generality in our representation of operators in divergence form, as it 
does not affect the result of Hill's method. 

Following the construction of the generalized periodic Evans function (13. 4j) . we may 
rewrite the truncated eigenvalue equation 

(3.9) {L a .j- XI)U = 
as 

(3.10) (I + Kj)U = 0, 

where K,j = JC^j + K, 2 ,j is the truncation of the Fourier representation K, = K,\ + JC2 of 
operator K to frequencies |(j, k)\ < J, that is, 

(3.11) = Vj(V 2 j - I)-U v and JC 2 ,j = (Vj - I)-\Ao,j + 1 - A). 

Continuing to follow the above construction of D a , we now define the truncated periodic 
Evans function as 

(3.12) D a>J (X) := det 2 (J - Kj) 
and notice that we have the following preliminary result. 

Lemma 3.7. The zeros of D a j correspond in location and multiplicity with those of L a j. 

Proof. This is immediate by the nonsingularity of (T>j — I)™ 1 and properties of the (usual, 
finite-dimensional) characteristic polynomial, together with the observation that 

deta(/ - Kj) = det 2 (Vj - I)- 1 det 2 (Vj + VjAi,j + A ,j - A/). 

□ 
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With this construction in hand, we now state the main result of this section. 

Theorem 3.8. For Aj G L^([0,X]), the sequence of determinants D a j converges to D a 
as J —> oo uniformly in a and locally uniformly in A. 

Proof. This convergence result follows from the proof of Theorem 13.41 Indeed, noting that 
D a> j is exactly such a sequence of approximate determinants, corresponding here to the 
ascending sequence of sinusoidal functions of integer wave number, by which the generalized 
periodic Evans function D a was defined in ()3.4f) . we find by our definition of the 2-modified 
Fredholm determinant that D a j — > D a pointwise in A as J — > oo for each fixed a G [0, 2tt). 
Moreover, recalling that the rate of convergence is determined by the difference between 
truncated operator tC j and IC in £>2 norm, and noting that we have uniformly bounded £>2 
estimates on each entry of /Cj, we find that this convergence is uniform in a and locally 
uniform in A. □ 

From Theorem 13.81 we immediately have convergence of Hill's method, as described in 
the introduction. For completeness, we state this result in the following corollary. 

Corollary 3.9. For Aj G L:L r ([0, X]), the eigenvalues of L a j defined in (|3.8|) approach the 
eigenvalues of L a in location and multiplicity as J — > oo, uniformly on |A| < R, a G [0, 2ir], 
for any R such that dB(0,R) contains no eigenvalues of L a . 

Proof. This is immediate from Theorem 13. 4^ Lemma [33 and Theorem 13. 81 along with basic 
properties of uniformly convergent analytic functions. □ 



3.3 Rates of Convergence 

Next, we address the issue of the rates of convergence of D a j to D a and of the approximate 
spectra to the exact spectra. Assuming slightly more regularity on the function A± in (|3.ip . 
we have the following easy convergence result. 

Theorem 3.10. For Aj G Hp eT ([0, X]) and each fixed R > 0, there exists a constant 
C = C(R) > such that for each fixed |A| < R 

\D <T , J (X)-D a (X)\<CJ- 1 / 2 . 

In particular, this estimate is locally uniform in A and uniform in a. 

Proof. The rate of convergence is bounded by \\KLj — IC\\s 2 from which we readily obtain 
the result using the Cauchy-Schwarz estimate 

E i^wi 2 ^ E i^'r 2 E ui 2 i^(i)i 2 < (c/j)\\A m \\ Hlmx]) 
m>j \j\>j \j\>j 

for each m G {0, 1}. For details, see the very similar estimates in the proof of Theorem 4.9, 
[GLZ] . □ 
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Notice that Theorem 13.101 does not imply a rate of convergence of the roots of D a j to 
the roots of D a , or, equivalently, the eigenvalues of L a j to the eigenvalues of L a . Indeed, 
the above convergence result is, with or without rate information, essentially an abstract 
one. For, though we find convergence the of analytic functions D a> j to D a , we don't obtain 
rates of convergence of their zeros without more structural information about D a itself. In 
particular, we can not conclude convergence rates of the approximate spectra to the true 
eigenvalues of L a using only the knowledge of the eigenvalues of L&,j computed in the course 
of Hill's method. This suggests the idea of computing the approximate Evans function D^j 
directly, instead of using it as a purely analytical tool, an idea that would be interesting for 
future investigation. Though in principle slower due to the need for multiple evaluations 
of eigenvalues, this computation is better conditioned, so there might perhaps be some 
counterbalancing advantages to this approach, besides the possibility already mentioned to 
obtain a posteriori estimates on the error bounds for eigenvalue approximations. We leave 
this as an interesting topic for further investigation, related to the larger question of relative 
advantages of standard periodic Evans function (as in [G] l vs. Hill's computations. 

4 Generalizations 

Here, we briefly discuss various generalizations of the theory developed in Section [3l 

4.1 Operators with nontrivial principal coefficient 

Consider now a system of the more general form 

(4.1) L a = dlA 2 + d x A 1 {a,x)+A {a,x), 

where A 2 is symmetric positive definite, satisfying A 2 {x) > C for some C > 0, uniformly 
on x 6 [0, X]. Define as usual A 2 to be the infinite-dimensional matrix representation of A 2 
under Fourier transform; that is, A 2 jk = A 2 (j — k). Then clearly A 2 is symmetric and, by 
Parseval's identity, satisfies A 2 > C when considered as a quadratic form on £ 2 (N). As a 
consequence, the J th truncation A 2) j, as a principal minor of a positive definite symmetric 
matrix, must also be positive definite and satisfy the same bound A 2 j > C. 
In particular, A 2 is invertible with 

A 2 X > 1/C, A^j > 1/C. 

Lemma 4.1. ||j4i3||g 2 < |A|£2 ||.B||,b 2; where \ ■ denotes L 2 ([0,X]) operator norm. 

Proof. Straightforward from the definition of || • ||b 2 . □ 

Corollary 4.2. For Aj S Lp Cr ([0,X]) and A 2 symmetric positive definite with A 2 (x) > C, 
the operator Ai := A 2 ^K, is Hilbert- Schmidt where K, = K,\+ K, 2 is defined as in (|3.1ip . 
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In this case, following the notation of Corollary 14.21 we define the generalized Evans 
function as D a (\) := det2(I — -A4), noting that the eigenvalue problem may be written 
equivalently as (I—M)U = 0. The associated series of Fredholm approximants is D CT j(X) := 
det2(I — Mj), with D CT) j(A) — > D a (\) uniformly as J — > oo, just as before, and zeros 
of D a corresponding in location and multiplicity with eigenvalues of L a . However, the 
corresponding object obtained by Hill's method is not the truncated Fredholm determinant 
T) (j.j defined above, but rather the modified version 

(4.2) D a>J (X) :=det 2 (I-^/Cj), 

and it is this function whose zeros correspond with the eigenvalues of the Hill approximant 
operator L a< j. 

To verify convergence of Hill's method in this case then, it is sufficient to show that 

(4.3) \\Mj - A^jKjWb, = \\{A^lC)j - A^jKj\\ B2 -> 

as J oo. Indeed, with this convergence result in hand we may conclude by f|2.3|) that 
limj^oo \D(j,J — D a j\ = 0, and thus D a j — > D a as J — > oo, yielding the convergence result 
as before. 

Theorem 4.3. For operators of the form (14. ip . Hill's method converges in location and 
multiplicity provided that Aj G Lp er ([0, X]). 

Proof. We sketch the proof of (|4.3p . By boundedness of H^IU^/rQ^n, we may truncate A2 
at wave number M to obtain an M-banded infinite-dimensional diagonal matrix centered 
around zero- frequency approximating A2 to arbitrarily small order in the ^ 2 (N) operator 
norm. Hence, for purposes of this argument, we may assume without loss of generality that 
A2 is M-banded diagonal operator centered about zero-frequency. Furthermore, noting 

that since A^ 1 is bounded in L 2 (R), for J G N sufficiently large the columns of A^ 1 
corresponding to frequencies \j\ < J — M are small off the principal 2 J + 1 — M minor and 
hence a brief calculation revealsthat 

(A 2 1 ) J A 2 ,j = 



where Em and Fm are M x M matrices that are invertible by invertibility of (A^^jAj, 
a property of principal minors of positive-definite symmetric matrices. By a further left- 
multiplication by the block-diagonal matrix 
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we obtain /2J+I) demonstrating that (A2,j)~ 1 agrees with (A~ 1 )2,j on the central 2J — 
2M + 1 dimensional minor. Recalling that \\JC — /Cj||b 2 — > as J — )• by (|3.5p . we thus 
obtain by a straightforward calculation 

11(^2 l K>)j ~ ^2J^j\\b 2 ~ 11(^2 ~ ^j^-^IIb 2 0) 

completing the proof by (|2.3j) □ 

4.2 Composite and Higher-order operators 

The reader may easily verify that all of the arguments of Sections [3] and 14.11 carry over to 
the case when the operator (jl.ip is replaced by a general periodic-coefficient operator 

L = d™a m (x) + d™- l a m -i{x) + ■■■ + a (x) 

where aj £ Lp er ([0, X]) and where the principal coefficient a m symmetric positive definite. 
Indeed, the analysis parallels that of previous sections except that one must substitute for 
{d 2 — 1) everywhere the positive definite symmetric Fourier multiplier 

\d 2 x - l\ m ' 2 = T~ l {\j\ 2 + l) m / 2 F, 

where j denotes the Fourier wave number and J- denotes Fourier transform. With these 
substitutions, our previous arguments immediately yield convergence of Hill's method in 
this case as well. 

Furthermore, it is straightforward to verify that all of the analysis of Sections 131 and I4TT1 
extends readily to the case of operators of "composite" type 

L = 

with a 3 k G Lp er ([0, X]) and a J mj symmetric positive definite for each suitable choice of indices: 
that is, still assuming L is a nondegenerate ordinary differential operator in some sense. 

Remark 4.4. It is the above observation that applies to the numerics in WJNRZ1\ {BJNRZ21. 
where the authors use Hill's method to numerically analyze the spectrum of the linearized 
St. Venant equations 

At - ct' - v! = 0, 

Xu - cv! - {f- z {F~ l - 2vu x )t)' = -{s + l)f s u r T - rf s+l u r ~ l u + v{f- 2 u')' 

about a given periodic or homoclinic orbit (u,f), where r, s, F, and v are physical param- 
eters in the problem and A is the corresponding spectral parameter. 
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4.3 Operators with general coefficients 

Our results are completely general in the scalar case, applying to all nondegenerate oper- 
ators. However, they are restricted in the system case by the condition that the principal 
coefficient(s) be symmetric positive definite. Whether this condition may be relaxed is an 
interesting operator-theoretic question regarding properties of Toeplitz matrices. 

Specifically, the property that we need to carry out Hill's method (and indeed, to 
complete our entire convergence analysis) is that the minor A.2 j of a Toeplitz matrix 
[^jmn = — n) be invertible for J sufficiently large. The question is what proper- 

ties of A2(x) are sufficient to guarantee this: in particular, is uniform invertibility enough? 
Alternatively, what are sufficient conditions on ^2? This seems an interesting problem for 
further investigation. 

Acknowledgement. Thanks to Bernard Deconink for pointing out the references 
[CuDllCDKKllDK] . 
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